Crossover from Intermittent to Continuum Dynamics for Locally Driven Colloids 
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We simulate a colloid with charge qd driven through a disordered assembly of interacting colloids 
with charge q and show that, for q^ ~ q, the velocity- force relation is nonlinear and the velocity 
fluctuations of the driven particle are highly intermittent with a 1// characteristic. When qd S> q, 
the average velocity drops, the velocity force relation becomes linear, and the velocity fluctuations 
are Gaussian. We discuss the results in terms of a crossover from strongly intermittent heterogeneous 
dynamics to continuum dynamics. We also make several predictions for the transient response in 
the different regimes. 
PACS numbers:82.70.Dd 



An individual particle driven through an overdamped 
medium exhibits a linear velocity vs applied force rela- 
tion. When quenched disorder is added to the system, a 
critical threshold driving force must be applied before the 
particle can move, and once motion begins, the velocity- 
force curves can be nonlinear. Examples of overdamped 
systems with quenched disorder that exhibit threshold 
forces and nonlinear velocity force curves include driven 
vortices in type-II superconductors or Josephson junc- 
tions [1], sliding charge density waves [2,3], driven mag- 
netic bubble arrays [4], and charge transport in assem- 
blies of metallic dots [5] . These systems have been stud- 
ied extensively and exhibit a variety of dynamic phases 
as well as power law velocity force curves. 

Another system that has been far less studied is an 
overdamped particle moving in the absence of quenched 
disorder but in the presence of a disordered two- 
component background of non-driven particles. Since 
there is no quenched disorder, simple overdamped mo- 
tion with an increased damping constant might be ex- 
pected. Instead, a critical threshold force F t h for motion 
exists and the velocity force curves are nonlinear with 
the power law form V — (F — Fth) 13 , as shown in recent 
simulations [6] and experiments [7] for individual colloids 
driven through a background of non-driven colloids. Un- 
like systems with quenched disorder, for F < F t h the 
entire system must move along with the driven particle. 
For F > F t h the driven particle is able to shear past its 
neighbors. The simulations give /3 = 1.5 for colloids in- 
teracting with a screened Coulomb potential. The experi- 
ments were performed with lightly charged colloids where 
steric interactions are important, and give » 1.5 in the 
lower density limit [7] . A particle moving through a vis- 
cous fluid can be regarded as interacting with many much 
smaller particles. In this limit, the surrounding medium 
can be replaced with a continuum and the velocity force 
curves are linear. In the case of a single driven colloid, it 
is not known when and how the dynamics change from 
nonlinear to linear, since it is expected that the system 
passes to the continuum limit when the colloids in the 
surrounding medium are small. This change could be a 
sharp transition, or it could occur as a cutoff of the scal- 



ing, or it could occur through a continuously changing 
exponent. 

In this work we consider a single colloid of charge qd 
driven through a T = disordered two-component as- 
sembly of other colloids with average charge q. We find 
that when q ks qd, the velocity force curves have a power 
law form with j3 > 1 that is robust over two decades 
and for different system sizes. In this regime the veloc- 
ity fluctuations of the driven colloid are highly intermit- 
tent and the colloid velocity V frequently drops nearly 
to zero when background colloids trap the driven colloid 
until rearrangements release it and V jumps back to a 
higher value. The fluctuations in V have a highly skewed 
distribution and 1// noise fluctuation properties. As qd 
increases for fixed drive Fd, the average velocity V drops, 
the velocity fluctuations become Gaussian, and /3 is re- 
duced. When qd 3> q, the driven colloid interacts with 
a large number of surrounding colloids and forms a cir- 
cular depletion region, while for qd <C q, the background 
colloids act as stationary disorder and the velocity force 
curves are linear. Our system can be experimentally re- 
alized for dielectric colloids driven with optical traps or 
magnetic colloids driven with external magnetic fields. 
In systems where it is difficult to vary the charge on 
individual particles, qd could be increased by capturing 
a large number of particles in a single optical trap and 
dragging the assembly through the background. Other 
related systems include dragging different sized particles 
through granular media [8] . There have been several pro- 
posals to use individual particle manipulation as a new 
microrheology method for examining frequency responses 
in soft matter systems [9]. It would be valuable to un- 
derstand under what conditions the individual particle is 
in the continuum overdamped regime or in the nonlinear 
regime in the dc limit. 

We consider a substrate-free, zero temperature, two- 
dimensional system with periodic boundary conditions 
in the x and y directions. A binary mixture of N = 
N c — 1 background colloids, charged with a ratio <?i/<72 = 
1/2, interact with a repulsive screened Yukawa potential, 
V{ r ij) = {li1j/\ r i ~ r J |)exp(-K|r l - r^), where q i{j) is 
the charge and r^-) is the position of colloid and 



1 



(a) (b) 
FIG. 1. Small dots: background of binary charged colloids. 
Large dot: Driven colloid with charge (a) q d /q = 1.33 and (b) 
qd/q = 67. 

is the screening length which is set to 2 in all our sim- 
ulations. Throughout the paper we refer to the average 
background charge q = (g x +^2) /2. The initial disordered 
configuration for the two-component background of col- 
loids is obtained by annealing from a high temperature. 
An additional driven colloid with charge qd is placed in 
the system and a constant driving force = F^x is ap- 
plied only to that colloid. The overdamped equation of 
motion for colloid i is 



<ir£ 
dt 



Ff + F d + F 7 



(1) 



where FJ 



nil 



V 



1, and the thermal 



force Ft comes from random Langevin kicks. We have 
considered various temperatures and discuss the T = 
case here. We have previously used similar Langevin dy- 
namics for colloids under noncquilibrium and equilibrium 
conditions [10]. The interaction range is assumed much 
larger than the physical particle size, and in this low vol- 
ume fraction limit, hydrodynamic interactions can be ne- 
glected and may be strongly screened [11]. To generate 
velocity-force curves, we set Fd to a fixed value and mea- 
sure the average velocity of the driven colloid < V > in 
the direction of the drive for several million time steps 
to ensure that a steady state is reached. The drive is 
then increased and the procedure repeated. Near the 
depinning threshold F d > F t h, the relative velocity fluc- 
tuations are strongly enhanced. In the absence of any 
other particles, the driven colloid moves at the velocity 
of the applied drive giving a linear velocity force curve. 
In this work we consider system sizes of L = 24, 36, and 
48 with a fixed colloid density of 1.1. This is four times 
denser than the system considered in Ref. [6]. 

In Fig. 1 we show images from the two limits of our sys- 
tem. In Fig. 1(a) the driven colloid (large black dot) has 
qd/q = 1-33 and is similar in charge to the colloids form- 
ing the surrounding disordered medium (smaller black 
dots). Fig. 1(b) illustrates the case qd/q = 67, where a 
large depletion zone forms around the driven colloid. 

To illustrate the scaling in the velocity force curves, we 
plot representative V vs Fd — F t h curves in Fig. 2 for 
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FIG. 2. Velocity V vs applied force Fd — Fth, where Fth 
is the threshold force for motion. From top to bottom, 
q d /q = 0.25, 1.33, 13, and 67. Dotted lines: fits with, from top 
to bottom, /3 = 1.0, 1.54, 1.28, and 1.13. System sizes: L = 24 
(open circles), 36 (plus signs), and 48 (filled diamonds). 

varied qd/q and different system sizes. Here F t h is the 
threshold velocity and the charge of the driven colloid 
increases from the top curve to the bottom. All of the 
curves have a power law velocity force scaling of the form 

V oc (F d - F th f (2) 

This scaling is robust over two decades in driving force. 
To test for finite size effects, we conducted simulations 
with systems of size L = 24, L = 36, and L = 48, indi- 
cated by different symbols in Fig. 2, and we find that the 
same scaling holds for all the system sizes. The scaling of 
the velocity force curves for the small charge q d /q — 0.25 
is linear, as seen in previous simulations [6]. In this 
regime the driven colloid does not cause any distortions 
in the surrounding media as it moves. This situation is 
very similar to a single particle moving in a quenched 
background, where it is known that the velocity force 
curves scale linearly or sublinearly [2]. As the charge of 
the driven colloid increases, the scaling exponent initially 
rises, as shown for the case of q d /q = 1.33 with (3 — 1.54, 
but the exponent decreases again for the more highly 
charged driven colloids, since q d /q — 13 gives (3 — 1.28 
and qd/q = 67 gives (3 — 1.13. As qd/q increases, the 
average velocity at fixed F d — F t h decreases when more 
background colloids become involved in the motion. 

We plot the changes in the scaling exponent (3 with 
varying q/q d from a series of simulations in Fig. 3, which 
shows three regions. For low qd/q, the motion is mainly 
elastic with (3 near 1. As q d /q approaches 1, the driven 
colloid charge becomes of the same order as that of the 
surrounding medium and the motion becomes plastic 
with (3 w 1.5. As q d /q increases further, (3 decreases 
approximately logarithmically toward 1, indicating that 
the motion of the surrounding medium is becoming more 
continuum-like. The maximum value of 
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FIG. 3. The scaling exponent (3 vs q d /q extracted from the 
velocity force curves. 

[3 falls at higher q d /q for less dense systems, such as that 
in Ref. [6], and at lower q d /q for denser systems. 

For driven colloids with qd/q -C 1, the background col- 
loids act as a stationary disorder potential. The driven 
colloid passes through this potential, deviating around 
background colloids as necessary, but the background col- 
loids do not respond to the presence of the driven par- 
ticle and remain essentially fixed in their locations. In 
contrast, when qd/q 1, the driven colloid strongly dis- 
torts the background and forms a depletion zone which 
moves with the driven colloid. Thus a large number of 
background colloids must rearrange in order to allow the 
driven colloid to pass, producing a continuum-like be- 
havior. Between these two limits, when qd/q ~ 1, no 
depletion zone forms, but when the driven colloid moves, 
it distorts the background which deforms plastically in 
order to allow the driven colloid to pass. Here we find in- 
termittent motion in which the driven colloid sometimes 
slips past a background colloid similar to the qd/q <C 1 
case, but at other becomes trapped behind a background 
colloid and pushes it over some distance, similar to the 
Qd/l^ 1 case. It is in the regime of this complex motion, 
when all of the charges are similar in magnitude, that the 
strongest deviation from linear response, [3 ~ 1.5, occurs. 

We next consider the velocity fluctuations of the driven 
colloid in the regime where (3 ~ 1.5 as well as in the high 
qd/q regime where (3 starts to approach 1. In Fig. 4(a) 
we plot a segment of the time series of the instantaneous 
driven colloid velocity for the case of qd/q — 1-33 at a 
drive producing an average velocity of V — 0.0425. The 
motion is highly intermittent and at times the colloid 
temporarily stops moving in the direction of the drive. 
When the driven colloid is trapped, strain accumulates in 
the surrounding media until one or more of the surround- 
ing colloids suddenly shifts by a distance larger than the 
average interparticle spacing and the driven particle be- 
gins to move again. As the drive is further increased, the 
length of the time intervals during which the driven 
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FIG. 4. (a) A segment of the instantaneous driven particle 
velocity time series V(t) for a system with q d /q = 1-33 at a 
drive giving an average velocity of V = 0.0425. (b) A segment 
of V(t) for qd/q = 67 for a drive giving the same average 
velocity as in (a), (c) The velocity distribution P(V) for the 
system shown in (a), (d) P(V) for the system in (b). (e) 
Power spectrum S(u) of the time series in (a). Dashed line: 
power law fit with exponent a — —0.8. (f) S(v) of the time 
series in (b). Dashed line: power law fit with a = —0.15. 

colloid is stopped decreases. 

In Fig. 4(b) we show V(t) for a system with qd/q = 67 
where the applied force gives the same average velocity as 
in Fig. 4(a). Here the amplitude of the velocity fluctua- 
tions is much smaller than the qd/q = 1-33 case and there 
are no intermittent stall periods. This strongly charged 
driven colloid is interacting with a much larger number 
of surrounding colloids than the weakly charged driven 
colloid would, and as a result, it cannot be trapped be- 
hind a single background colloid, giving much smoother 
motion. We note that we find no intermittent behavior 
for the strongly charged driven colloid even at the low- 
est applied forces. We also measured the variance of the 
transverse velocity fluctuations V yi and find that it de- 
creases with increasing qd/q roughly as a power law with 
an exponent of -1.7. This decrease occurs since a larger 
number of background colloids are contributing to the 
fluctuations experienced by the driven particle, leading 
to a smoother signal. 

In Fig. 4(c) we plot the histogram of the velocity fluc- 
tuations P(V) for the time series shown in Fig. 4(a). The 
fluctuations are non-Gaussian and are heavily skewed to- 
ward the positive velocities with a spike at V = due 
to the intermittency. We note that in simulations of vor- 
tex systems where nonlinear velocity force scaling occurs, 
bimodal velocity distributions are also observed when in- 
dividual vortices are intermittently pinned for a period 
of time before moving again [12]. Non-Gaussian velocity 
fluctuations are also found in sheared dusty plasmas [13]. 
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FIG. 5. The transient velocity response to a sudden applied 
force of F d /Fth = 0.8. (a) q d /q = 1.33; the dashed line is a 
power law fit with an exponent of a = —1.1. (b) q d /q — 67; 
the dashed line is an exponential fit. 



For higher drives, we find that the average velocity in- 
creases; however, the histograms remain highly skewed 
for all charge and drive regimes where the scaling in the 
velocity force curves give a large [3 ~ 1.5. For compar- 
ison, in Fig. 4(d) we plot P(V) for the large q d /q sys- 
tem shown in Fig. 4(b). Here the histogram has very 
little skewness and fits well to a Gaussian distribution. 
For other drives at this charge ratio we observe similar 
Gaussian distributions of the velocity. In general, we find 
decreasing skewness in the velocity distributions as qd/q 
increases. The Gaussian nature of the fluctuations is also 
consistent with the interpretation that, as q d /q becomes 
large, the system enters the continuum limit. 

In Fig. 4(e), we show that the power spectrum S{v) for 
the time series in Fig. 4(a) has al//" form with a = 0.8. 
Throughout the (3 « 1.5 regime we find similar spectra 
with a = 0.5 to 1.1, indicative of intermittent dynamics. 
For comparison, in Fig. 4(f) we show that S(v) for the 
high qd/q case has a white velocity spectrum character- 
istic, indicative of the absence of long time correlations 
in the velocity. In general, a decreases with increasing 
q d /q. For the small charge regime of qd/q <C 1, where the 
driven colloid moves without distorting the background, 
we also observe a white noise spectrum. Here the velocity 
fluctuations are determined by the static configuration of 
the background particles, and P(V) does not show any 
intermittent periods of zero velocity because if the col- 
loid becomes trapped in this regime, there can be no 
rearrangements of the surrounding medium to untrap it. 

To explore the transient behavior of the system, we 
consider the effect of a suddenly applied subthreshold 
drive of Fd/F t h — 0.8. In Fig. 5(a) we show the transient 
velocity response for a system with qd/q — 1.33. Here the 
velocity relaxation is consistent with a power law decay, 
V(t) oc t~ 1A . The driven colloid translates by several 
lattice constants before coming to rest with respect to 
the surrounding medium. The large velocity oscillations 
that appear at longer times are due to the local plastic 
rearrangements of the surrounding colloids as the driven 
colloid passes. We find a power law decay in the transient 
response for values of q d /q that give (3 > 1.28. If the 
suddenly applied subthreshold drive is small enough that 



no local rearrangements of the surrounding medium are 
possible, then an exponential decay of the velocity occurs 
instead. In Fig. 5(b) we show the transient response for 
the case of qd/q = 67. Here the relaxation is fit to V(t) ex 
exp(— t). The large velocity fluctuations that appeared 
for the smaller q d are absent. For all the large values of 
qd/q we observe an exponential velocity relaxation, and 
we also find exponential relaxation for very small charges 
qd/q<- 1.0. 

In summary, we have studied a single colloid with vary- 
ing charge driven through a disordered background of 
other colloids in the absence of quenched disorder. When 
the charge of the driven colloid is close to the same 
as that of the surrounding colloids, a nonlinear power 
law velocity force curve appears with an exponent near 
(3 = 1.5. In this regime, the time dependent velocity 
is intermittent with a highly skewed velocity distribu- 
tion and l/f a noise fluctuations. As the charge of the 
driven colloid is increased, the velocity force characteris- 
tic becomes more linear while the effective damping from 
the background colloids increases. The number of col- 
loids that interact with the driven colloid increases and 
the velocity fluctuations become Gaussian with a white 
noise spectrum. We predict that in the nonlinear regime, 
the transient velocity responses are of a power law form, 
while in the linear regime the transient velocity is expo- 
nentially damped. We interpret our results as a crossover 
in the response of the background colloids from intermit- 
tent dynamics when the driven and background colloids 
are similarly charged to continuum dynamics for a highly 
charged driven colloid. 
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